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Abstract 

We study a fermionic dark matter in a non-supersymmetric extension of the standard model 
with a family symmetry based on Dq x Z2 x Z2. In our model, the final state of the dark 
matter annihilation is determined to be e~^e~ by the flavor symmetry, which is consistent with 
the PAMELA result. At flrst, we show that our dark matter mass should be within the range 
of 230 GeV — 750 GeV in the WMAP analysis combined with — )• 67 constraint. Moreover 
we simultaneously explain the experiments of direct and indirect detection, by simply adding 
a gauge and D^ singlet real scalar fleld. In the direct detection experiments, we show that the 
lighter dark matter mass ~ 230 GeV and the lighter standard model Higgs boson ~ 115 GeV is 
in favor of the observed bounds reported by CDMS II and XENONIOO. In the indirect detection 
experiments, we explain the positron excess reported by PAMELA through the Breit-Wigner 
enhancement mechanism. We also show that our model is consistent with no antiproton excess 
suggested by PAMELA. 
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1 Introduction 



The existence of the dark matter (DM) in the Universe has been estabhshed by measurements. The 
WMAP experiment tells us that the amount of the DM is considered about 23% of energy density 
of the Universe [T]. As indirect detection experiments of the DM, PAMELA reported excess of 
positron fraction in the cosmic ray [2]. This observation can be explained by annihilation and/or 
decay of DM particles with mass of O{10'^'^) GeV. The PAMELA experiment searches antiproton 
as well in the cosmic ray, and it is consistent with the background [3] . Therefore, if these signals are 
from annihilation and/or decay processes of DM particles, this indicates that the leptophilic DM is 
preferable. However, even if the DM is leptophilic, the resultant positron fraction depends on the 
generation of final state leptons. For instance, if the final state of annihilation and/or decay of the 
DM is r+r", it will overproduce gamma-rays as final state radiation Therefore it is considerable 
that the leptophilic DM can refiect fiavor structure of elementary particles. In this point of view, 
several works discussing the DM and fiavor structure have been done so far [5| 1^ 171 [5| 1^ [TU | [TT | [T ^ [T3] . 

Flavor structure of elementary particles is thought to be determined by symmetry, so called fiavor 
symmetry [l^. In our previous work [13], we have discussed fermionic DM model with the standard 
model (SM) extension with the Dq fiavor symmetry [TB]. In this model, three generations of matter 
fields including right-handed neutrinos are embedded into doublet and singlet representations of Dq 
group in particular way. The light neutrino masses are induced by radiative correction through 
inert SU{2)l doublet Higgs bosons r] which do not have vacuum expectation values(VEVs) [T6| [T7]. 
We identify a heavy Majorana neutrino of Dq singlet ns with the DM candidate. The DM is 
stable because of the additional Z2 symmetry. Since the Dq symmetry completely determines fiavor 
structure of the model, the final states of annihilation of the DM via Dirac Yukawa interaction Vj'^Lns 
are fixed to be electron-positron pair and r neutrino pair. In that paper, we have found that the DM 
mass is constrained to be in the range 230 GeV — 750 GeV from the condition of the relic abundance 
and — )■ 67 process. However, this annihilation of the DM via ?7-mediated t- and u-channel processes 
does not give enough s-wave contribution to the cross section because it is proportional to mass of 
the final state irte^v Therefore, this model requires very large enhancement (~ 10^) of the cross 
section at the present Universe compared with that at the early Universe to explain the PAMELA 
data, which is not realistic. 

In this paper, we extend our Dq model of Ref . |13j by adding gauge and Dq singlet scalar field 
(p, which couples with the DM ns as (fusus- While the final states of the DM annihilation are the 
same as those of the previous model, there exist s-channel annihilation processes mediated by (f. 
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In this case, the Breit-Wigner enhancement mechanism works which can give enough boost factor 
P^[T^ I2U]. Moreover, the new field mixes with the singlet Higgs doublet (psi which is responsible 
for mass of the quark sector. This mixing can simultaneously induce antiproton production by DM 
annihilation and interaction with quarks in atoms. We find that the spin-independent cross section 
of the DM and quarks via the mixing between ip and (ps can be close to sensitivities of direct DM 
detection experiments such as CDMS II [21j and XENONIOO [22], suppressing antiproton fiux in the 
cosmic ray. 

This paper is organized as follows. In section 2, we review our model briefiy and summarize 
the predictions for lepton sector coming from the fiavor symmetry. In section 3, we analyze the 
Higgs potential and mixing between the SM Higgs and new singlet scalar ip. In section 4, we show 
constraints of DM mass from WMAP and /i — > 67 process. We discuss direct and indirect detection 
of DM in section 5 and 6, respectively. Section 7 is devoted to conclusions and discussions. 

2 The Model 

In this section, we briefiy review a SM extension with x Z2 x Z2 family symmetry [13]. 
2.1 Yukawa couplings 

We introduce three "generations" of Higgs doublets 4>i^s, inert doublets r^/ 5, and one generation of 
inert singlet (p. Where / = 1,2 and 5* denote Dq doublet and singlet, respectively, and assume that 
each field is charged in specific way under the family symmetry shown in Table [1] and [21 
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Table 1: The Dq x Z2 x Z2 assignment for the leptons. The subscript S indicates a Dq singlet, and 
the subscript / running from 1 to 2 stands for a Dq doublet. Lj and Ls denote the S'f/( 2) /^-doublet 
leptons, while ej, e|, rij and ns are the 5'f/(2)i-singlet leptons. 
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Table 2: The Dq x Z2 x Z2 assignment for the Higgs bosons. 



Under the Z2 symmetry (which plays the role of R parity in the MSSM), only the right-handed 
neutrinos nj,ns and the inert Higgs doublets rji^rjs are odd. All quarks are assumed to be singlet 
under the family symmetry so that the quark sector is basically the same as the SM, where the Dq 
singlet Higgs doublet (ps with (+, +) of Z2 x Z2 plays a role in the SM Higgs in the quark sector. No 
other Higgs bosons can couple to the quark sector at the tree-level. In this way we can avoid tree- 
level flavor changing neutral currents (FCNCs) in the quark sector. The Z2 symmetry is introduced 
to forbid tree-level couplings of the Dq singlet Higgs (ps with L/, L5, rii and n^, simultaneously to 
forbid tree- level couplings of 0/, 77/ and rjs with quarks. As shall be discussed later, the gauge singlet 
if plays an important role in explaining an indirect detection reported by PAMELA. Furthermore, 
it is expected to explain the direct detection as CDMS II, because our dark matter rig, singlet 
right-handed neutrino, couples to the quark sector by small mixing between Lp and 0s, which should 
be estimated to satisfy the experimental results. We will show the numerical analysis of the mixing 
for both experiments later. 

The most general renormalizable DqX Z2X. Z2 invariant Yukawa interactions in the lepton sector 
are found to be 

C-Y = Yl [Y:,\L,ta2<P,)el + Y:MLa)nt] (2.1) 

a,b,d=l,2,S 

- 2^ —niui Y^sns - —^pniuj - —LpusUs + h.c, (2.2) 

1=1,2 1=1,2 

where the coupling constants ©i,s are complex in general. The electroweak symmetry is broken by 
the VEVs (0i) = (02) = vn/2 , (0s) = vs/V2, V = vl + vl = (246 GeV)^ (r/,,s) = (v^) = 
and we obtain the following mass matrix Me and diagonalization matrix Uf.L of MgM]^ in the 



3 



charged lepton sector: 
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where m| 4 5 are real parameters whose values are determined by observed charged lepton masses 

and = m^/m,-. In the neutrino 
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sector, Yukawa couplings in the mass eigenstates are given by 
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where the Dirac Yukawa couplings hi {i = 2,3,4) are of order one. Notice that the Dq singlet 
right-handed neutrino couples only with Ls and rjs- Since we consider the case that rjj^s are 
inert doublets which do not have VEVs, Dirac neutrino mass matrix is not generated and canoni- 
cal seesaw mechanism does not work. Light Majorana neutrino masses are generated by radiative 
seesaw mechanism at one-loop level [H]. In this mechanism, Majorana mass is proportional to 
hf kV'^ M / {M"^ — m"^)), where k denotes typical coupling constant of non self-adjoint terms in 
the Higgs potential. When k = 0, an exact lepton number U{1)l' invariance is recovered, where the 
right-handed neutrinos n/ ^ are neutral under U{1)l' in contrast to the conventional seesaw models. 
This U{1)l/ forbids the neutrino masses, so that the smallness of the neutrino masses has a natural 
meaning. Now we can derive some predictions of our model based on the family symmetry: 

1. If ee,^ = 0, the mixing matrix UeL has the maximal mixing in its right-upper block which is the 
origin of the maximal mixing of atmospheric neutrino mixing. Only an inverted mass spectrum 
rui^s < f^fi 2 is allowed. 



2. Non- zero 6*13 is predicted as sin^ 9 
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1.2 X 10 ^. This small value of 61^ is consistent with 



the best fit value 0.020+|]:[|J]^ with la error 



3. The effective Majorana mass (m)ee is bounded from below as {m)ee ~ 0.02 eV. 

As a result of this discussion, we can assume that Mi 5 = 0{TeV), k and hi = 0{1). Moreover, 
as can been seen from Eq. fl2.5p . if one identifies the Dq singlet right-handed neutrino ns to be the 
DM, it mainly couples with electron (and positron) with large coupling /13 ~ 1. This selection rule is 
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remarkably determined by the family symmetry. These facts play a crucial role in the study of cold 
DM (CDM) as discussed below. 



3 Higgs Potential 



In this section, we analyze the Higgs potential. As discussed in Refs.[T3l [23], the Higgs potential 
consists of Dq symmetric and breaking terms. Since Dq invariant Higgs potential has an accidental 
global 0(2) symmetry, the latter must be introduced in order to forbid massless Nambu-Goldstone 
(NG) bosons. Essentially, such soft breaking terms are mass terms of the Higgs bosons. For the 
potential of (0/, (ps), the soft Dq breaking mass terms [23] are given by 



+ 



+ h.c. 



(3.1) 



where is real while ^\ is complex in general. The mass term of (0/, ips) is dominated by Eq. fl3.ip . 
and subdominantly given by invariant terms of order V'^. One finds that the breaking terms 
Eq. (l3.ip preserve the minimum symmetry 5*2 under 0i ^ 02- The key point is that the 5*2 invariance 
is required not only to ensure the vacuum alignment (0i) = (02) 7^ but also to forbid NG bosons 
which violate the electroweak precision test of the SM. 

Since the Higgs potential of 0/^5 and rji^s are analyzed in Ref.[T3], we do not explicitly show that 
here again. In the present model, the new field is introduced and it plays an important role in our 
analysis. Therefore we explicitly show the potential including ip. The most general renormalizable 
Dq X Z2 X Z2 invariant Higgs potential of (f is given by 

V{ip) = mlip + mlip^ + m3(p^ + Xiip'^, (3.2) 

V{(f),(p) = m4(0t.0s)v5 + m5(0|07)v9 + A2(0^05)V5^ + A3(0|0/)V3^ (3.3) 

Viv,^) = V{<P,^)i<P^v), (3.4) 

where all parameters are considered to be real without loss of generality. By using the decomposition 
of SU{2)l doublets 0/^5, 



Pi + ^o-/ 
V2 I 7207 

we find the mass matrix of neutral Higgs bosons as 
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where p = [pi, ps), o" = (ct/, as). Each 3x3 element M^^^'s are given by [13] 
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where the coefficients of 0(1). The (/^-dependent terms are given by 



The stable minimum conditions are found by partially differentiating the potential by (p as 
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Therefore, we obtain the vacuum conditions for {(pi^s) ^ind {(p) = as 



m 



1 + o i^s^^ + ^0^5) = 0, MJ^ > 0, vs(D)m4(5) > 



(3.14) 



The mass matrix M| is diagonalized by the 7x7 orthogonal matrix O, as OMlO"^ . Notice that 
quarks couple only with 05 via Yukawa interactions, and (psips) mixes with (f via m^. This mixing 
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parameter will induce both interaction of the DM with atoms (direct detection) and antiproton 
flux in the cosmic ray. We will discuss these DM phenomenology below. Note also that there is no 
mixing between and r] because 777,5 do not get VEVs. 

The SM Higgs is described in terms of the linear combination of flavor eigenstate fields as 

SM - Higgs = Onpi + C12P2 + Onps + O^a^ + dscis + Oieas + Oircp, (3.15) 

and the other combinations correspond to heavy neutral Higgs bosons with mass of several hundred 
GeV. Therefore the ps — f mixing is proportional to Oj^Oj^. In the following analysis, we give 
numerical values of the matrix O. 



4 WMAP and /i ^ 67 Constraint 

In this section, we derive conditions for mass of the DM Ms and charged component of t] boson M, 
following the result of Ref . [13] . 



4.1 fi ^ e-j Constraint 

The DM mass Ms is constrained from the fi ^ ej process. The branching fraction of yU — )■ 67 from 
Fig{T]is approximately given by 

r2 /T\/r2\ 



B(M - n) = J^A- ^ |X^900 GeV^I^ , A- . ,4 (4.1) 

and 



^, , 1 — 6x + Sx^ + — 6a;2 In a; , 

^^<^' = EiT^- • <^-2> 

where x = M^/M^. As we can see from the Yukawa matrix of Eq.f l2.5p . only r]s couples to ns with ei 
and with ^l, where the coupling with /i^ is suppressed by rrie/m^ ~ 0.005. In the next subsection, 
we will obtain the constraints of the DM mass Ms which is consistent with the observed DM relic 
density fi^/i^ ~ 0.12 pQ and fi — 67, assuming to be the DM. 



4.2 WMAP 

In the analysis of Ref. |13j. we have found that it is more natural and promising that only ns of 
three right-handed neutrinos remains as a fermionic CDM candidate. Furthermore since charged 
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Figure 1: The contribution to the — )■ 67 process. 



component of rjs boson couples to cl and ns due to our original matrix in Eq. (l2.5p . it remarkably 
leads to be a clean signal if the charged extra Higgs boson rjs is produced at LHC. 

We simply find the thermally averaged cross section {(Jiv) for the annihilation of two n^'s [25] 
from Figj2] in the limit of the vanishing final state lepton masses: 



fli + bi 
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X 



ai = 0, bi 



Ms 



/i3|V(l -2r + 2H 
247rM2 



(4.3) 
(4.4) 



where is 775 mass, Ms isns mass which is our DM candidate and T is temperature of the Universe. 
The thermally averaged cross section Eg. (14. 3 p does not contain s-wave contribution as a consequence 



ns 



I 
I 



+ 



ns 




Figure 2: Annihilation diagrams of ns for the cross section aiv. 

of massless limit of the final state particles, and we find that the allowed region for the DM mass is 
around (9(10^) GeV from the constraints of WMAP results [1] and — )■ 67 decay. 

In FigJ3]we present the allowed region in the — Ms plane, in which Vt^h? = 0.12 and — )• 
67) < 1.2 X 10~^^ [26] are satisfied, where we take {h^l < 1.5. From Eg. (14.30 . retaining = 0{1) 
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is quite important to find the promising DM mass regions, as we mentioned before. Note tliat tliere 
is no allowed region even for {h^l ^ 0.8. As can been seen from FigJHl we find tlie mass range as 
follows: 

230 GeV < Ms < 750 GeV, 300 GeV < < 750 GeV. (4.5) 

In this analysis, we have calculated the mass bound for Sunyaev and Zeldovich (SZ) effect [27]. In our 
model, 7]g, which decays to high energy e^, may affect the CMB by the inverse Compton scattering, 
if the lifetime is not between 10"^^"^-* sec. From the condition that the lifetime of rjg comes into the 
allowed region, mass has the bound of 30 GeV < < 750 GeV. Where the Yukawa coupling 
nearly equals to 1, and ^ Ms are assumed. Hence, one finds that the SZ effect satisfies the both 
constraints of /i — )■ 67 and cosmological pair annihilation of CDMs sufficiently. 
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Figure 3: The allowed region in the - Ms plane in which fi^/i^ = 0.12, 67) < 1.2 x 10"^^ and {h^l < 1.5 

are satisfied. 

5 Direct Detection 

We analyze the direct detection search through the experiments of CDMS II [21] and XENONIOO 
[22] . The main contribution to the spin-independent cross section is from the t-channel diagram with 
the mixing between ip and (ps, as depicted in FigJH Then the resultant cross section for a proton is 
given by 
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Figure 4: The t-channel diagram for the direct detection [21] by the if — (ps mixing. 



with the hadronic matrix element 

A = y /(f)^ + l y (5.2) 

^ q=u,d,s ^ q=c,b,t ^ 

where rup is the proton mass. The effective vertex ag in our case is given by 

<^,^^Sf(^. ,,3, 

'"'SM-Higgs 

Here rusM-Higgs is the SM Higgs mass and V oc {nig/v) is a Yukawa couphng constant of the quark 
sector. Notice that the quark sector couples only to (ps- In the numerical analysis, we set the Higgs 
masses to avoid the lepton flavor violation (LFV) process as follows: 

115 GeV < rnsM-Higgs < 200 GeV, 500 GeV < other six neutral Higgs boson masses. (5.4) 

Under this setup, the elastic cross section is shown in FigJSl Where we set lO^^^Of^Sgl = 0.1, which 
we call the "mixing". We plot the DM mass Ms in the region 10 — 1000 GeV. Since the allowed 
region of the DM mass is 230 GeV — 750 GeV from the WMAP analysis combined with /i — )■ 67 
constraint, rather smaller SM Higgs mass is favored if these experiments could detect the DM near 
the current bound. 



6 Indirect Detection 

The PAMELA experiment implies that there could be positron excess [2], but not be antiproton 
excess [3j. In order to describe the PAMELA results successfully through an annihilation process of 
the DM, we need enhancement of the cross section by using the Breit-Wigner mechanism [TH] . 
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Figure 5: The spin-independent cross section as a function of the DM mass for the direct detection 
[211 [22] • The "mixing" is defined by \O^^Oj^es\ and set to be 0.1. The longitudinal black line 
represents the SM Higgs boson mass range. 



6.1 Positron Production from DM annihilation 

The main channel of the DM annihilation in the present Universe is depicted in FigEl The ns 
annihilation cross section to leptons is given by 
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(6.1) 

(6.2) 
(6.3) 
(6.4) 
(6.5) 
(6.6) 



where the spin of initial states is averag ed and a = M^/M^ and /3 = ml/M^. Notice that Eq.(EI]) 
has the s-wave contributions because the coupling ©s is complex . The mass parameter is a mass 
eigenvalue of the Higgs mass matrix which is satisfied the resonance mass relation Mr ~ 2mJ3, 



"'^We take account of physical pole (A > 0). Unphysical pole analysis is studied in detail in [TB] 
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Figure 6: The main process for the positron excess from the DM annihilation. The s-channel diagram induces the 
Brcit-Wigner enhancement. 



Tji is the decay width to r;,5r;,5G and A = 1 — 4M|/M^. The resonance particle R is described in 
terms of the hnear combination of flavor eigenstate fields as 

R = OriPi + Or2P2 + Or^ps + Or^cji + C^5^2 + ORf^as + Or^^. (6.7) 

There are the other contributions to the ns annihilation cross section such as t, w-channel in Figj2]or 
the interference contributions between t, u-channel and s-channel. However all we have to consider 
is the contribution of Eq.f l6.ip because this is dominant at the present Universe that the DM relative 
velocity is f ~ 10~'^. One finds that the flavor symmetry remarkably fixes the final states to be 
positron/electron in our scenario 

The thermally averaged annihilation cross section (cr2f ) is defined as 

(^zt') = r-PT—r^ — TiFTio — ) (6-8) 



f d^Pid^P2frf'^ 



eg 
2 



where pi is the momentum of initial particle i and = e~^^^'^ is the Maxwell-Boltzmann distribution 
function. If we can expand the annihilation cross section in terms of f ^ as = 02 + ^2'^^, we can 
calculate it easily as (o"2f) = 02 + 662/a^, where x = Ms/T. Although such a naive treatment is 
not justified when the annihilation cross section has a resonance point, an approximate estimation 
is obtained as follows if the condition 7/j ^ A is satisfied [HI [20] : 

1 ^2 .fX 

x^/^e--^, (6.9) 



^ ^ {ihlHORrflm&sY ^ ^2 , , ^ ^2 .2I ^A 



(47r)% 



4„l/2 



(Re6,)M^ + |)+(Im6.)^/| 



K 7« 



^Although there are other decay channels like 4's4>S: we assume that decay to nsns is dominant to lead to the 
Breit-Wigner enhancement. 

•^We assume that nj in the loop does not contribute to the positron production because nj can produce the tauon 
final state with no suppression, which is now forbidden by the Fermi-LAT 7-ray experiment [i]. Such a condition can 
be realized in our model by controlling the coupling 61 to be small. 
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Figure 7: Contours of the boost factor (BF) in the 6s complex plane. The parameter A is taken as 10 in the left 
figure and 10^^'' in the right figure. We set the other parameters as Mjj — 500 GeV, Ms = 230GeV and Ihsl'^Oj^'^ — 1. 



where 'jn = Th/Mr. Since is proportional to v A, one might suspect that large annihilation cross 
section is obtained under the condition 7ij/A <^ 1. We will discuss this point below. 
We define the boost factor BF as 

BF^ (6-10) 

3.0 X 10-9[GeV~Y V ; 

and contours of the boost factor are shown in FigJTl where the red regions satisfy the condition 
•Jr/ A < 0.1 and = 500 GeV is taken as a typical example. The degree of the fine tuning is 
smaller (larger) if the smaller (larger) value is taken because the thermally averaged cross section 
{a'2v) is inversely proportional to M^. One finds that a large boost factor is obtained through the 
Breit-Wigner enhancement from FigjT] if the parameters satisfy A < 10^^'^ and Im&s ^ ReSs. 
Under this condition, the thermally averaged annihilation cross section and the decay width of R are 
written as 

16v^, 



(An) 



Ihl'Oj,, (ReGs) 



2 2 



7712 
ri 



.3/2^ 



-xA 



iOTT 



(6.11) 



(6.12) 



thus one find that ImSg ^ Re©5 is important to obtain large annihilation cross section and small 
decay width. 

The fiux of positron and electron from DM annihilation is given by $e± (e) [2S] and the positron 
fraction is given by 



Positron Fraction 



$e+(e) + $!f -(e) + $e-(e) + $ 



prim. 



e) + 



(6.13) 
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where (e) are the contributions from DM annihilation, and the others are the background fluxes 
given by 

^pnm./ X u^^ut fGeV^^cm^^s^^sr^M 

$«^^-fe) = rii^^^ fGeV^^cm-^s-^sr"^) 

$r+-(e) = j^fGeV-^cm-^s-^sr-i). (6.14) 




Figure 8: The positron fraction for DM annihilation(ns) into e+e . The red, green and blue hnes are the best fits 
for Ms = 230 GeV and (craw) = 8.5 x IQ-^ GeY'^, Ms = 450 GeV and (ctzv) = 2.6 x IQ-^ GeV"^ and Ms = 750 GeV 
and {(J2v) ~ 6.8 x 10^^ GeV~^ respectively. 

The direct positron fraction is plotted in Figj8]for some fixed parameters. The BF of order 10^ 
is required in all cases. This BF is not large enough to fit the Fermi-LAT data [29j. Thus the 
constraints from diffuse gamma rays and neutrinos are not severe as long as isothermal dark matter 
profile is considered [4j. It might be worth mentioning that the DM mass less than O (TeV) is in 
favor of the experiment recently reported by HESS [30], if one considers the NFW profile [3T| . 

6.2 Muon Flux Measurement from Super-Kamiokande 

We briefly mention that the high energy neutrinos induced by DM annihilations in the earth, the 
sun, and the galactic center are an important signal for the indirect detection of the DM [32]. Such 
energetic neutrinos induce upward through-going muons from charged current interactions, which 
provide the most effective signatures in Super-Kamiokande (SK) [33j. Once the thermally averaged 
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cross section of the muon flux reaches the same order of the cross section required by the PAMELA 
results, it is natural to expect that such a value of cross section is close to the upper bound of the 
muon flux measured by SK. In fact, our model has the large cross section enhanced by the Breit- 
Wigner mechanism with Vj. pair final state as can been seen in Figj6l However since the total cross 
section is proportional to the neutrino mass as in Eq.f l6.ip . the neutrino flux is extremely suppressed 
than positron and electron fluxes. 

6.3 Antiproton Production from DM Annihilation 




0.1 1 10 100 

T [GeV] 



Figure 9: The antiproton flux is mainly induced by the top-pair production through the Higgs mixing between Lp 
and 4>S in our model. Where T = E — Ms is the kinetic energy of antiproton, we set {a^v) = 3 x 10~^ GeV~^. 



Finally, we briefly discuss the antiproton flux in the cosmic ray. Since our model has the quark- 
DM coupling through the Higgs mixing between ip and 0s, as discussed in section [5l we have to verify 
that our antiproton flux is consistent with the antiproton experiment of PAMELA. The main source 
comes from the top quark pair production, and substantially bottom and charm pair production. 
The cross section of the nsris — )• qq processes is given by 



(^■iv{,nsns ^ qq) = — Wl 

/TT 



Ml vl {s - Mlf + MlTl 



R 



X 



{OmO 



R7) 



Ml 



{OmO 



R7) 



(6.15) 



where the index q is summed over top, bottom, and charm quark. The energy-squared s of the initial 
state is defined in Eq. (16.51) . 
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The thermally averaged annihilation cross section is expressed in terms of A, 75 and some cou- 
plings as 



{Im&sY I ml ml \ ( ml 



167ri/2 V MlvlMl 



The PAMELA experiment implies the positron excess, but no antiproton excess. Thus the ratio of 
the annihilation cross section to leptons and quarks constrains the mixing parameters between <f and 
(ps- The ratio is given by 



M \mj V ^1^1/ 1^3 P (Re65) + (111165) 

where /i, I2 and /a are taken as (9(10~^) which is evaluated by numerical analysis. If we require the 
boost factors for leptons and quarks to be 100 and 1 respectively, the constraint to the couplings 
becomes 

^^l ^ ^""^ 2 ^ 0{lQ-^'), (6.18) 

(Re65)' + (Im©^)' 

where we have taken the masses of Ms = 450 GeV and = 500 GeV. We find that the mixing 
matrix elements Ors and Orq which appear in {asv) need to be suppressed by C(10~^^) in order to 
have no antiproton excess if ©5 is 0{1). 

The fiux of antiproton from DM annihilation is given in Ref . ^H] • We plot the antiproton flux as a 
function of the kinetic energy of antiproton T = E — Ms in Figj9j Where we adopt {(J3V) = 3 x 10~^ 
GeV~^, i.e. BF = 1, which is required to explain the WMAP experiment, and the same set up as 
the positron case. The key parameters contributing to the direct detection are O13 and Ou which 
come from SM Higgs mediation, while those to the indirect detection of the antiproton are 0_r3, 
Orq, and Ofn which come from resonant bosons. It suggests that both of them can be explained by 
independent way. Hence it is easy to find the allowed region avoiding such an enhancement as well 
by controlling many parameters in the Higgs sector. 



7 Summary and Conclusions 

In this paper, we have considered that two important issues of the dark matter in a non-supersymmetric 
extension of the radiative seesaw model with a family symmetry based on Dq x Z2 x Z2: direct de- 
tection recently reported by CDMS II and indirect detection reported by PAMELA. We suppose 
that the Dq singlet right-handed neutrino is the promising candidate of the DM. Analyzing the 
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H ^ e'f together with the WMAP result, we have shown the allowed region for the DM mass to 
be 230 GeV < Ms < 750 GeV, within a perturbative regime. In the analysis of the direct de- 
tection experiment of CDMS II and XENONIOO, we have shown that the Higgs mixing between 
ip and (ps plays an important role in generating the quark effective couplings, and also there exist 
allowed region to be detected by those experiments in near future. As a result of the positron pro- 
duction analysis through PAMELA, a couple of remarks are in order. In the case of Ms = 230 GeV, 
Ms = 450 GeV and Ms = 750 GeV, each of {av) = 8.5 x 10-^ GeV-^ {av) = 2.6 x 10"^ GeV"^ and 
(cTf ) = 6.8 X 10~^ GeV^^ is required, respectively. In all cases the required boost factor is at most 
~ 0(10^), which is realized by the Breit-Wigner enhancement mechanism if the sutable parameter 
regions are choosen. Also such boost factor is not large enough to fit the Fermi-LAT result. Thus 
constraints from diffuse gamma rays and neutrinos are not severe as long as isothermal dark matter 
profile is considered. Finally, we have investigated the antiproton flux in the cosmic ray to compare 
to the direct detection. We found that the constraint of the mixing from the direct detection can 
easily satisfy the allowed region for no antiproton excess by controlling many parameters in the Higgs 
sector. 
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